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ABSTRACT 

Optical/UV radiation from accretion disks in quasars is likely to be partly scattered by a hot plasma enveloping 
the disk. We investigate whether the scattering may produce the steep rises in polarization observed blueward of 
the Lyman limit in some quasars. We suggest and assess two models. In the first model, primary disk radiation 
with a Lyman edge in absorption passes through a static ionized "skin" covering the disk, which has a temperature 
kT ~ 3 keV and a Thomson optical depth tt ~ 1. Electron scattering in the skin smears out the edge and produces 
a steep rise in polarization at A < 912 A. In the second model, the scattering occurs in a hot coronal plasma 
outflowing from the disk with a mildly relativistic velocity. We find that the second model better explains the data. 
The ability of the models to fit the observed rises in polarization is illustrated with the quasar PG 1630+377. 

Subject headings: accretion, accretion disks - galaxies: active - polarization - quasars: individual 
(PG 1630+377) - radiative transfer - ultraviolet: galaxies 
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1. INTRODUCTION 

Steep rises in polarization blueward of the Lyman limit have 
been detected in 3 of 10 observed quasars (Impey et al. 1995; 
Koratkar et al. 1995). Especially challenging is the case of 
PG 1630+377 where the polarization increases from ~ 2% at 
1000 A to — 20% at 700 A. Attempts to explain such a steep 
rise have not been successful (e.g., Blaes & Agol 1996; see 
Koratkar & Blaes 1999 for a review). One toy model fitting 
the data was suggested by Shields, Wobus, & Husfeld (1998). 
They assume, however, an arbitrary large ad hoc jump in polar- 
ization at A = 912 A in the disk comoving frame and then fit 
the observed rise by the relativistically smeared-out jump. 

A possible physical mechanism producing the rises in po- 
larization is Compton upscattering of the disk radiation by a 
hot plasma. Assuming that the disk emits UV radiation with a 
Lyman edge in absorption, one finds that the upscattering (1) 
smears out the Lyman edge and (2) polarizes the radiation. The 
upscattered radiation dominates the spectrum blueward of the 
edge, and this leads to the rise in polarization at A < 912 A. 

We investigate here two scenarios: 

A. Multiple scattering in a static slab with Thomson optical 
depth tt ~ 1. The scattered radiation is limb-brightened in the 
slab and the multiply scattered radiation acquires a polarization 
parallel to the disk normal (Sunyaev & Titarchuk 1985). 

B. Scattering in a mildly relativistic outflow of optical depth 
tt < 1. Relativistic aberration causes a limb-brightening of 
the disk radiation in the outflow rest frame. This effect leads 
to a parallel polarization, p ~ 20%, of the scattered radiation 
(Beloborodov 1998, hereafter B98). 

Transfer of polarized radiation in a static slab corona of a 
high temperature, kT ~ 100 keV, was previously studied in de- 
tail (e.g., Haardt & Matt 1993; Poutanen & Svensson 1996; Hsu 
& Blaes 1998). It has been shown that the corona produces par- 
allel polarization rising toward short wavelengths. Possible rel- 
evance to the UV polarization rises in QSOs has been assessed 
by Hsu & Blaes (1998) who concluded that the predicted rise is 
not steep enough to explain the data. 
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In this Letter, we argue that the scattering in a layer at a 
Compton temperature, kTc ~ 3 keV, would generate a much 
steeper rise in polarization as compared to a hot static corona 
with kT ~ 100 keV. The layer may be associated with an ion- 
ized "skin" of an accretion disk heated by the corona emission. 
In §2, we perform calculations of radiative transfer in the skin 
demonstrating the production of steep rises in polarization. 

Then, in §3, we investigate the outflow model. Previously, it 
was suggested as a possible mechanism producing parallel op- 
tical polarization (B98). We now find that the model can repro- 
duce the rises in UV polarization. The results are summarized 
in §4. 

2. COMPTONIZATION IN A STATIC SLAB 

2. 1 . Production of parallel polarization 

Let semi-isotropic radiation propagate into an ionized slab 
from an underlying disk. With increasing height, the attenu- 
ated intensity of unscattered radiation, 1$ (the 0-th scattering 
order), gets suppressed at large angles with respect to the nor- 
mal because the slab has large optical depths at these angles. 
The limb-darkening of Iq results in perpendicular polarization 
in the first scattering order. In the second order, photons come 
from the slab itself, and their intensity, I±, is limb-brightened 
if the slab optical depth tt £ I. The radiation scattered twice, 
I2, then acquires a parallel polarization. The limb-brightening 
progresses with successive scatterings and then saturates as the 
intensity Jjv approaches the eigen-function of the transfer prob- 
lem (Sunyaev & Titarchuk 1985). The resulting parallel polar- 
ization progressively increases with scattering orders and then 
saturates. For tt ~ 1, the saturation occurs after ~ 5 scatter- 
ings. 

In a hot slab, the scattering is accompanied by a photon 
blueshift. After N scatterings, a photon of initial energy e ac- 
quires an energy of- (l + 4kT/m e c 2 ) N e (see, e.g., Rybicki & 
Lightman 1979). One thus expects a strong impact of the slab 
on both the disk spectrum and the polarization. If the original 
UV spectrum has a Lyman edge in absorption then the edge 
energy, £l = 13.6 eV, delineates two regions in the emerging 
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spectrum. At e > el, the upscattered radiation dominates the 
observed flux, and the further away from the edge, the higher 
are the scattering orders contributing to the spectrum. The ra- 
diation is therefore parallelly polarized at e » e L . By contrast, 
at £ < El, all scattering orders contribute to the spectrum, and 
this leads to a small net polarization. 

A slab with temperature kT <~ 100 keV produces a grad- 
ual growth in parallel polarization toward high energies, so 
that p reaches - 10 - 15% at e > 1 keV (Hsu & Blaes 
1998). The polarization maximum is blueshifted far away from 
the edge because the gain in energy per scattering is large, 
As/e <~ AkT/m e c 2 — 4/5. The high T is thus the reason 
for the slow growth of p with respect to wavelength. If the 
Comptonization occurred in a slab of a relatively low temper- 
ature, then a large number of scatterings would be accompa- 
nied by a modest blueshift, and p(A) at A < 912 A would 
be much steeper. With this motivation, we consider a transi- 
tion layer ("skin") between the UV disk and the corona as a 
possible location of the Comptonization. The layer is Comp- 
ton heated by the corona emission to an equilibrium tempera- 
ture kTc of a few keV and has an optical depth tt ~ 1 (e.g., 
Zycki et al. 1994). The Kompaneets' y-parameter of the layer 
is 4(fcT '/ 'm e c 2 )T^ <C 1, which implies a very soft spectrum of 
upscattered photons. 

If the disk has a patchy hot corona then the bulk of UV ra- 
diation comes directly from the disk without being scattered 
in the corona (Haardt, Maraschi, & Ghisellini 1994). At the 
same time, the whole UV disk may be covered by the ionized 
skin. Then the skin may dominate the Comptonized spectrum 
at modest energies, while the patchy corona dominates the hard 
X-ray spectrum. The skin may contribute to the soft excess ob- 
served in quasar spectra. 

We model the skin as a slab of optical depth tt ~ 1 and 
temperature kT ~ 1 — 10 keV. Below the slab, we assume a 
blackbody source with temperature kT hh — 3 eV. The Lyman 
edge in the source spectrum is modeled as a sharp jump down 
to zero intensity at 912 A. In the calculations of radiative trans- 
fer in the slab, we use the iterative scattering method described 
in detail in Poutanen & Svensson (1996). 

2.2. Results 

The results for tt = 1 and kT =1, 3, and 10 keV are pre- 
sented in Figure 1 for a disk inclination cos# = 0.23. The 
polarization is close to its maximum at this inclination. With 
increasing tt, p decreases, changing its sign, and in the limit 
tt 3> 1 one arrives at the standard Chandrasekhar-Sobolev 
perpendicular polarization. An increase in T leads to a larger 
energy gain per scattering that makes the polarization rise less 
steep. On the other hand, a decrease in T and/or tt leads to 
a reduction of the flux blueward of the Lyman edge. This is a 
general tendency of a Comptonization model: dramatic rises in 
polarization imply modest fluxes blueward of the edge. A com- 
promise is achieved when tt ~ 1 and kT ~ 3 keV. In this case, 
the predicted behavior of the flux and polarization resembles 
the data on PG 1630+377 (see Fig. 1). The edge is smeared-out 
at large inclinations, while it is strongly pronounced at mod- 
est inclinations (see Fig. 2). Additional smearing-out may be 
present due to relativistic effects near a Kerr black hole. 

A robust feature of the model is that p changes sign across the 
edge, i.e., the small polarization redward of the edge is perpen- 
dicular to the disk normal (see also Haardt & Matt 1993). By 
contrast, the data favor the same position angle blueward and 
redward of the edge (see Koratkar et al. 1995). The constancy 



of the position angle is better seen in PG 1222+228 (Impey et al. 
1995). Besides, optical polarization measured in many objects 
tends to be parallel to the disk normal (e.g., Antonucci 1992). 
The predicted perpendicular redward polarization thus appears 
to contradict the data. 
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FIG. 1 . — Intensity, I\ (in arbitrary units), and polarization, p, from a disk 
covered by an ionized "skin" of optical depth r T = 1, as viewed at incli- 
nation 9 = 11°. The sign of p is chosen so that p > corresponds to the 
case where the electric vector is parallel to the disk normal. The dashed, solid, 
and long-dashed curves correspond to skin temperatures kT = 1, 3, and 10 
keV, respectively. Below the skin, we assume unpolarized blackbody emission 
of temperature 3 eV with a sharp Lyman edge in absorption. The data for PG 
1630+377 (Koratkar et al. 1995) show the detections of p above the 2.5<r level. 
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FIG. 2. — Skin model with tt = 1 and kT = 3 keV. Top: Emerging 
intensity at different inclinations, p, = cos 0. The dashed, solid, long-dashed, 
and dotted curves correspond to p, = 0.047, 0.23, 0.5, and 0.95, respectively. 
Bottom: Angular dependence of polarization at different wavelengths. 
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3. SCATTERING IN A CORONAL OUTFLOW 

A corona above an accretion disk need not be static and for- 
mation of mildly relativistic outflows is very likely. Moreover, 
there are indications from X-ray observations for bulk motion in 
coronas of black hole sources (Beloborodov 1999a; Zdziarski, 
Lubihski, & Smith 1999). Here, we consider a toy model in 
which the outflow is replaced by a slab of outflowing plasma 
with a given column density, N c . The underlying UV disk is 
assumed to emit isotropically. A characteristic velocity of the 
outflow is f3 = v/c ~ /3* = (4 - V7)/3 w 0.45, which cor- 
responds to equilibrium with the disk radiation field (i.e., the 
bulk acceleration by the radiative pressure is balanced by the 
Compton drag at /? = /?*, Gurevich & Rumyantsev 1965). The 
equilibrium must be established in an outflow composed of 
pairs (Beloborodov 1999b). When viewed from the outflow co- 
moving frame, the disk radiation is limb-brightened. This effect 
has a maximum at (3 = /?* (B98). 

A typical polarization of the scattered radiation is ~ 10 — 
20% depending on the disk inclination, /i (B98). The observed 
weak polarization redward of the Lyman edge, p 2%, indi- 
cates that the contribution of the scattered radiation to the ob- 
served flux is small. In a slab geometry, it implies that the out- 
flow is optically thin, N c < 1 . At the same time, blueward of 
the edge, the scattered radiation can dominate the observed flux, 
producing a steep rise in polarization at A < 912 A. If the pri- 
mary radiation crucially decreases at A < 912 A then the flux 
here is mainly composed of scattered photons that have been 
blueshifted upon scattering due to Doppler effect. In addition 
to the Doppler effect associated with bulk motion, there may 
be Doppler shifts due to thermal motions in the outflow. We 
therefore consider outflows with a range of temperatures, T. 

3.1. Calculations 

The equations of Thomson transfer of a polarized frequency- 
integrated radiation in a cold outflow are given in B98. Here, 
we extend the problem: (1) We include the dependence on fre- 
quency, v, and study the polarization along with the radiation 
spectrum. (2) We relax the assumption kT <C m e c 2 and calcu- 
late radiative transfer in outflows of high temperatures. 

In a slab geometry, the polarized radiation is described by the 
two Stokes parameters, and Q(v,fM). The scattering is 

represented by the source functions, Si(v, /i) and Sq(v, (j). We 
hereafter use vector notation I = (I, Q) and S = (Si, Sq). To 
calculate S, we transform I into the outflow comoving frame. 
Polarization, p = Q/I, is invariant with respect to Lorentz 
boosts along the normal, and I transforms as (see, e.g., Rybicki 
& Lightman 1979) 

l c {v c ,^ c ) = D s (n)l(v,ii), 



H-0 



v c = 7(1 - f3jj)v = D{n)v, 



1-/V 

where the index c stands for the comoving frame and 7 = 
(1 — /3 2 )~ 1//2 . In the comoving frame, we compute the source 
functions using the 2 x 2 redistribution matrix, R, derived for 
a hot static plasma (Nagirner & Poutanen 1993; see also Pouta- 
nen & Svensson 1996), 



We calculate the observed radiation in the single-scattering 
approximation, i.e., I b s ~ I + r M (S — I), where r p = 
N c <jt (1 — (3fJ,)/n is the slab optical depth at inclination fi. This 
approximation is reasonable when (1) the outflow has a small 
optical depth (77 < 1) and (2) we are interested in the radiation 
near the Lyman edge. The high scattering orders will affect 
the results at energies above a typical energy of once-scattered 
photons. We assume isotropic unpolarized primary radiation, I, 
(the primary depolarization in the disk is likely to be caused by 
Faraday rotation, e.g., Agol & Blaes 1996). 

3.2. Results 

First, we consider the case of a sharp Lyman edge (Fig. 3). 
After passing through the outflow, a polarized tail appears in the 
spectrum blueward of the edge. One can see from Figure 3 how 
a high outflow temperature increases the flux at A < 912 A. On 
the other hand, a high T reduces the polarization of the upscat- 
tered radiation (see also Poutanen & Vilhu 1993; Nagirner & 
Poutanen 1994). This effect is due to mildly relativistic thermal 
motions that cause random aberration of light in the electron 
rest frame. The random character of the "thermal aberration" 
leads to a reduction of the systematic polarization. 
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FIG. 3. — Disk radiation after passing through an outflow with /3 = /?* = 
0.45 and column density N c = 0.2ov^ , as viewed at fi = f3. I\ = Iv 1 /c 
is plotted in arbitrary units. The sign of p is chosen so that p > corresponds 
to the parallel polarization. The primary (unpolarized, isotropic) radiation of 
the disk is assumed to have a blackbody spectrum, fcT^b = 3 eV, with a sharp 
Lyman edge (dotted curve). The dashed, solid, and long-dashed curves repre- 
sent the observed radiation (in the single-scattering approximation) for outflow 
temperatures kT = 0, 50, and 100 keV, respectively. 



The amplitude of the tail is proportional to the outflow col- 
umn density. The discontinuity in the observed spectrum comes 
from the edge in the primary radiation. The polarization is only 
due to scattered radiation which has no discontinuity, hence 
Qobs is smooth across the edge. This gives a relation between 



and then transform S c into the lab frame, S(v, ji) 
L>- 3 0u)S c (z/ c ,Mc). 



bs: Pobs/Pobs = W-Cbs- Where +/- 

correspond to A just above/below 912 A. 

In order to emulate the observed spectrum of PG 1630+377, 
we use a toy primary intensity I £ oc e 2 cxp[— (e/£ ) 6 ] and 
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choose £o = 13 eV. This spectrum emulates a smeared-out Ly- 
man edge (shown by the dotted curve in Fig. 4). The original 
smearing-out may be the result of relativistic effects near the 
black hole. The emerging radiation after passing through the 
outflow is shown in Figure 4 for N c = 0.23ovp 1 at inclination 
6 = 60°. A reasonably good fit to the PG 1630+377 data is 
obtained for kT w 50 - 100 keV. 
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FIG. 4. — Outflow model with /3 = /?* = 0.45 and column density 
N c = 0.230'^, as viewed at inclination 6 = 60°, against the data on PG 
1630+377. The assumed original disk emission is shown by the dotted curve. 
The dashed, solid, and long-dashed curves correspond to outflow temperatures 
kT = 3, 50, and 100 keV, respectively. 

A special feature of the outflow model is that the maximum 
p is predicted at [i = (i which corresponds to [i c = 0, albeit the 
polarization is also high for a wide range of inclinations (see 
B98). More importantly, the parallel orientation is predicted at 
all wavelengths, in contrast to our first (static) model. Note also 
that the outflow model provides a steeper rise in p and predicts 



a hardening in the spectrum at ~ 750 A similar to that present 
in the data. 

4. SUMMARY 

We have found that Comptonization of UV radiation with a 
Lyman edge in absorption can produce very steep polarization 
rises blueward of the edge. Redward, the predicted polarization 
is low and weakly dependent on wavelength. 

We model the Comptonizing plasma as a hot slab with three 
parameters: (1) the optical depth, tt, (2) the temperature, T, 
and (3) the bulk velocity, [3. Our first (static) model emulates 
an ionized skin of an accretion disk. It has tt ~ 1, kT <~ 3 keV, 
and /3 = 0. The skin produces a steep rise in parallel polariza- 
tion blueward of the Lyman edge. The predicted polarization, 
however, changes its sign across the edge, which is not favored 
by the observations. Another potential difficulty is the strong 
Lyman edge predicted at modest inclinations. 

In our second (outflow) model, the hot plasma is optically 
thin, tt < 1, and its bulk velocity, j3 ~ 0.5, is in equilibrium 
with the disk radiation field (B98). A steep rise in polarization 
is then generated for a wide range of parameters. The predicted 
polarization is parallel to the disk normal at all wavelengths. 
Note that the redward polarization may be further reduced if an 
additional unpolarized source contributes at large wavelengths 
(e.g., the disk emission from larger radii, which is not scattered 
in the central outflow). 

Both models predict polarization rises only in those objects 
where there is a Lyman edge or a dramatic spectral break in the 
original disk spectrum. A large polarization, p <; 20%, is then 
predicted at favorable inclinations. The skin model has maxi- 
mum p at cos 9 s=s 0.25, and the outflow model - at cos = (3. 

The outflow model appears to be more promising. Note that 
the slab geometry is a rough approximation for an outflow. A 
real outflow is at least 2D (axisymmetric) or it consists of lo- 
calized ejecta from the disk (Beloborodov 1999a). A full 2D 
model should include self-consistently the effects of relativis- 
tic disk rotation and gravitational redshift, which are important 
near a Kerr black hole. 

This work was supported by the Swedish Natural Science 
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